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m ' Abstract 

New solutions for static non-rotating thin disks of finite radius with 
nonzero radial stress are studied. A method to introduce either radial 
pressure or radial tension is presented. The method is based on the 
use of conformal transformations. 
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1 Introduction 



Axially symmetric solutions of Einstein field equations corresponding to disk- 
like configurations of matter are of great astrophysical interest and have been 
extensively studied. These solutions can be static or stationary and with or 
without radial pressure. Solutions for static disks without radial pressure 
were first studied by Bonnor and Sackfield |T|, and Morgan and Morgan [^], 
and with radial pressure by Morgan and Morgan ||, and, in connection with 
gravitational collapse, by Chamorro, Gregory and Stewart f|. Static disks 
with radial tension may also be of interest and, although their relevance in 
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astrophysics is limited, these solutions can be studied as possible sources of 
exact solutions for axially symmetric static problems. 

Static thin disks without radial stress generate a Weyl spacetime. Then 
they are described by two different metric functions |5], The stability 
of these models can be explained by either assuming the existence of hoop 
stresses or that the particles on the disk plane move under the action of their 
own gravitational field in such a way that as many particles move clockwise 
as counterclockwise. This last interpretation is frequently made since it can 
be invoked to mimic true rotational effects. A large class of statics thin 
disks solutions were obtained by Letelier and Oliveira using the inverse 
scattering method. Solutions for self-similar static disks were analyzed by 
Lynden-Bell and Pineault M, and Lemos M. The superposition of static 



disks with black holes were considered by Lemos and Letelier [10, 11, 12 



and Klein |T3[. Recently Bicak, Lynden-Bell and Katz JT-1 studied static 



disks as sources of known vacuum spacetimes and Bicak, Lynden-Bell and 
Pichon [31| found an infinity number of new static solutions. 

Static disks with radial stress are described by three metric functions so 
presenting some degree of complexity and, despite its clear relevance, has 
been seldom treated in the literature. The fact that the additional metric 
function can be considered as the real part of an analytic function will be the 
main ingredient of our method to introduce radial stress. Static disks with 
radial stress can also be obtained in other ways. Lemos || studied self-similar 
static disks working in spherical polar coordinates and using self-similarity 
to reduce the Einstein field equations to an ordinary differential equation. 
Recently, Pichon and Lynden-Bell |16|] obtained static disks by making a cut 
through a Weyl solution above all singularities or sources. The identification 
of this solution with its mirror image yields disks with radial stress. 

In this work we consider static non-rotating thin disks of finite radius 
with nonzero radial stress. To obtain a well behaved radial stress we find 
an analytical function compatible with the boundary conditions imposed 
on the additional metric function by the requirement that the solution be 
singularity-free outside the disk. We find a general expression for the radial 
stress that can be written in terms of a square root in such away that by 
choosing a sign solutions with radial pressure and solutions with radial ten- 
sion can be obtained. When we take the positive (negative) sign of the square 
root, the stress is positive (negative) and so we will have disks solutions with 
radial pressure (tension). 



2 



The plan of the paper is as follows. In Sec. 2 we present the principal 
aspects of static thin disks solutions with radial stress, the expression for the 
energy-momentum tensor of the disk and the general solution that leads to 
disks with radial pressure or with radial tension. In Sec. 3 we study solutions 
with radial pressure. In Sec. 4 solutions with radial tension are considered. 
We conclude, in Sec. 5, by summaring our main results. 



2 Static Thin Disks Solutions 
2.1 The General Solution 

The metric for a static axially symmetric spacetime can be written in quasi- 
cylindrical coordinates (r, ip, z) in the form |T7| , TB \ 

ds > = - e * d f + i?V* dtp 2 + e A ~*(dr 2 + dz 2 ), (1) 

where R, $, and A are functions of r and z only. In the vacuum, the Einstein 
equations are equivalent to 

R, rr + R, zz = 0, (2) 

(R$,r), r +(R$,z),z=0, (3) 

(R 2 +R 2 )A, r = \RR )T ($,2 -$,2 ) + RR, Z $, r + (R* +Ri ), r , (4) 

(R 2 +R 2 )A, Z = \RR 1Z ($,2 -$, r 2 ) + RR, r $, r + (i?, r 2 +R* ), z , (5) 

where ( ), a = d/dx a . 

In order to obtain a solution of @ - (|5|) representing a thin disk located 
at z = 0, < r < 1, we assume that the metric functions R, $, and A are 
continuous across the disk, but have discontinuous first derivatives in the 
direction normal to the disk. Althoug the radius of the disk has been taken 
as unity, a suitable re-scaling of r can be made. We will comeback to this 
point later. The reflectional symmetry of (|2|) - (||) with respect to the plane 
z = allows us to assume that R, $, and A are even functions of z. Hence, 
R, z , $, 2 and A, z are odd functions of z. We shall require that they not vanish 
on the surfaces z = ± . 
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We now solve the Laplace equation 
the boundary conditions 



R{0,z) 



0. 



in the region r > 0, z > with 



z > 0, 



(6) 



R, z {r,0) 



r > 1, 



(7) 



R, z (r,0) ^ 0, r < 1. 



(8) 



Let = r + iz. We can consider i? as the real part of an analytical function 
W{uj) = R(r, z) + iZ{r, z). Due to the discontinuous behavior of R, z across 
the disk, the Cauchy-Riemann equations imply that Z is discontinuous across 
the disk while its normal derivative Z, z is continuous. Then, the invariance 
of reflection in the plane of the disk allows us to assume W(r, —z) = W(r, z). 
We also require that W{uj) be locally one-to-one in the outer region of the 
disk. 

A simple solution compatible with these conditions is 



W(u) 



CO 



(9) 



where a is a positive constant, and the ± means that we can take any of the 
two possible signs of the square root. When z > we obtain 



R(r, z) 



Z(r, z) 



r ± a 



z ± a 



(|r 2 - z 2 - 1\ 2 + 4rV)§ + (r 2 - z 2 - V 



(\r 2 - z 2 - 1| 2 + 4r 2 z 2 )2 - (r 2 - z 



(10) 



11' 



For z < we only need to replace z — > — z in (|I0|) and (|Tl|). 

In order to solve the other field equations, we make a conformal transfor- 
mation 

r^R(r,z), z^Z{r,z), (12) 

in such a way that 



$ = $, A^A = A- l n \W 



i\2 



(13) 
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where W = dW/duj. With this transformation the field equations take the 
usual Weyl form |5|, [J 

®,rr + + ®,zz= 0, (14) 

A[$] = lJR [($, 2 R -$, 2 z )dR + 2$, fl ^dZ]. (15) 

Note that (|14|) is the Laplace equation in flat three-dimensional space and 
so $ can be taken as a solution of Laplace equation for an appropriated 
Newtonian source with axial symmetry. Once a solution <3> is known, A is 
computed from (fl5|) and so we obtain, from flL?!) - (fl3|), a solution of the field 
equations @ - ©■ 

By using ( Jl3| ) - (|15|) we find the nonzero components of the Riemann 
curvature tensor, 

Rrrnp = ?^[R,r& 2 z - $ ,\) + 2R, Z ® , R ® , z ], (16) 
= ^^[i?, z (l>,| -$,^)-2 J R, r l>, R <l, z ], (17) 

Rzipztp Rrtpztpi ( 1 ^) 

and the Kretschmann scalar 

i^e^ = x [Jtf - <£,| ) 2 

+ +6$,|$,| +$,|) (19) 

+ 4R, r R, g $, R $, z &l - $,|)]. 

From these expressions, using (P) and the properties of $, we can test if the 
solutions have strong curvature singularities. 

2.2 The Energy-Momentum Tensor 

Due to the discontinuous behavior of the derivatives of the metric tensor 
across the disk, the Riemann curvature tensor contain Dirac delta functions. 
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The energy-momentum tensor can be obtained by the distributional approach 
due to Papapetrou and Hamouni Hl9| , Lichnerowicz [p0|1 , and Taub See 
also Israel |Z5|, P5| . It can be written as T a b = [T a b ] S(z), where S is the 
Dirac function with support on the disk and [T a b ] is the distributional energy- 
momentum tensor, which yield the volume energy density and the principal 
stresses. The "true" surface energy-momentum tensor of the disk can be 
written as r a b = e^'^' 2 [T a b ]. 

The discontinuities in the first derivatives of the metric tensor can be cast 



as 



11, 24 



bab 9ab,z\ z=Q+ 9ab,z\ z=Q -i 

then the distributional energy-momentum tensor is given by 



(20) 



[T\] = -{b az 6\ - b zz 6\ + g az b z b - g zz b\ + b c c (g zz 6\ - g az 6 z b )}. 

For the metric ([!]) we obtain [T a b ] = diag(—e, p v , p r , 0), with 

2R, Z ~ 



Py 



2$„-A„ 



[A„ 



R 



(21) 

(22) 
(23) 



Pr 



2R, 



R 



(24) 



where all these quantities are evaluated at z = + , < r < 1. In order to 
obtain expressions for a disk with non-unit radius, we only need to make the 
transformation r — > ar, where a is the radius of the disk. 

From ([24]) and (§) we can compute the general expression for the radial 

stress p r , 

2ay/l - r 2 , 



Pr 



± 



i + {a 2 - iy 



(25) 



where the ± correspond to the two signs of the square root in (|9]). Thus, 
when we take the positive (negative) sign the stress is positive (negative), 
we have disks solutions with radial pressure (tension). From the fact that $ 
and A are both real and finite at the surface of the shells, we can see that p r 
is a well behaved function of r for any value of a. 
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By doing z = 0, < r < lin (|1(]) and ( |TTD we obtain the image of the 
disks by the conformal mapping, 

a 2 R 2 + Z 2 = a 2 , (26) 

so that the disks are mapped into spheroidal thin shells of matter. According 
with this, the solution leads to three different boundary conditions for a 
three-dimensional potential problem. These boundary conditions correspond 
to the three possibilities for a: a = 1, a spherical shell, a > 1, a prolate 
spheroidal shell, and < a < 1, an oblate spheroidal shell. 

Using the above shells as sources to solve the Laplace equation ( |14D we 
obtain many different solutions of (0) - (0). In all the cases the solutions 
of (|14|) are known and can be written in terms of Legendre functions. Due 
to the discontinuous character of the Weyl coordinate Z across the disk, a 
polynomial in odd powers of Z is discontinuous across the disk while that 
a polynomial in even powers is a continuous function everywhere but with 
discontinuous first derivatives with respect to z across the disk. Accordingly, 



the solutions of fll4| ) must be chosen with the appropriated Legendre functions 
of even order in such away that the solutions of (Q) - (|5|) have the assumed 
behavior at the disks. 

In order to have a one-to-one correspondence in the outer region to the 
disks, we need to choose the appropriate sign of the square root in (|9|) that 
map this region completely into the exterior or interior region to the shells. 
The appropriate choose will depend of the value of a. When a ^ 1 we only 
can take the positive sign of the square root because with the negative sign 
the map (§) is not single valued at the exterior of the disk and, in this case, 
we only can take the solutions of ([14]) that correspond to the potential in the 
exterior of the appropriated spheroidal shell. 

On the other hand, when a = 1 we can take separately both of the 
possibilities in (Q) and, for each one, we will obtain a well posed boundary 
value problem. When we take the positive sign, the outer region to the disk 
is mapped into the outer region to a spherical shell of radius one. Thus 
we must take the solutions of flUD that correspond to the potential in the 
exterior of the shell. When we take the negative sign, the outer region to the 



disk is mapped into the interior of the same shell and so the solutions of (|14]) 
must be those that correspond to the interior potential. 
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3 Static Disks with Radial Pressure 



For the positive sign of the square root in (^|) there are three different families 
of solutions corresponding to the three possibilities a = 1, a < 1 and a > 
1. In each case we consider simple known Weyl solutions and from them, 
through ( |12"D and (p~3|) , we obtain the corresponding thin disks solutions. 
As we shall see, in all the solutions considered we obtain disks with a well 
behaved central sector where the energy-momentum tensor agree with the 
weak and strong energy conditions P3 , but with a border region where e < 



in violation of the weak energy condition. 
3.1 Chazy-Curzon Disks 

The first family of solutions is obtained by choosing a = 1 and transforming 
the disk into a spherical shell of radius one and the outer region of the disk 
into the outer region to the shell. To solve the Laplace equation ( UM we 
introduce spherical coordinates (11,9), related to the Weyl coordinates by 

R = Ksm9, Z = 1Zcos9, (27) 

with < K < oo and -f < 9 < f . The shell is at 1Z = 1. 



The general solution for the Laplace equation (|T4|) in the exterior of the 
shell, continuous across the disk, can be written as 

- ~ 2C 2n P 2n (cos9) 

$ - - 2_ — — > i 28 ) 

n=0 /V 

where C 2n are constants and P 2n (cos9) are the usual Legendre polynomials. 
The simplest case is the Chazy-Curzon solution KB], 



s 2 7 7 7 2 sm 2 6' , . 

® = -^r, A = — — , 29 

n iz 2 ' v ; 

where 7 is a real constant. Note that $ and A are well behaved continuous 
functions for IZ > 1 and vanish at infinity. 

The metric tensor in Weyl coordinates is given by 

9,1 = -'Atwtw}' m 



8 



g w = Rexpl \ , (31) 




^{7f^-W?w)' (32) 

9zz = 9rr- (33) 



Using (10), (|TT]) and (13), we can obtain its expression in the original co- 
ordinates. Then, taking z = + and using (|22|) - (|24D, (0) and ([□]), we 
obtain 

e = p r [2(7-1) -7V], (34) 

Pv = Pr [1+ 7 V], (35) 

p r = p Vl — r 2 exp{7 2 r 2 }, (36) 

where po is a positive constant and < r < 1. 

In order to analyze the behavior of the above solution we consider the 
weak and strong energy conditions P5| . One of the requirements of the 
weak energy condition is that e > everywhere. From fl3"l|) we see that this 
condition is equivalent to 

2( 7 - 1) - 7 V > 0, 0<r<l. (37) 
This equation has one zero at 



2(7-1) 



V = ^-r^ (38) 

and so, in order to have e > for r = 0, we must impose the condition 7 > 1. 
The maximum value for ro is obtained for 7 = 2 and is given by tq 2 = 1/2. 

The strong energy condition require that cr > 0, where a = e + p v + p r is 
the effective Newtonian density. From fl3H|) - (|55|) we see that cr = 27p r and 
so is positive everywhere on the disk for all value of 7 > 0. In Figure 1 we 
plot the functions a, e, p r and p^ in the interval < r < 1 for this choice 
of the parameters, in po units. The disk have a central region that satisfy 
the weak and strong energy conditions and a border region where e < 0, in 
violation of the weak energy condition. 
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3.2 Zipoy-Voorhees Disks 

The second family of solutions is obtained by choosing a > 1 in (Q), trans- 
forming the disk into a prolate spheroidal shell. To solve the Laplace equation 
( |HD we introduce prolate spheroidal coordinates (£,??), related to the Weyl 
coordinates by 



R 2 



(39) 



where n = \fa 2 — 1, 1 < £ < oo and < r\ < 1. The shell is located at 
£ = a/Va 2 - 1 > 1- 

The general solution for the exterior of the shell, continuous across the 
disk, can be written as 



$ = -£2C 2n P 2n (r?)Q 2n (0, 



(40) 



n=0 



where Q 2n are the Legendre functions of second kind. The simplest case is 
the Zipoy-Voorhees solution [|28], p9"| , also known as the Weyl 7-solution |], 
given by 



$ = 7 In 



£ + 1' 



A = Y In ■ 



(41) 



£2_ ^2' 

where 7 is a real constant. We see that $ and A are continuous on the disk 
and vanish at infinity. 

The metric tensor in the prolate spheroidal coordinates is given by 



9tt 



R 2 



£ - 1 



£ + 1 



9z 



\W 



9r 



>\2 





7 2 




£, 2 ~ v 2 







(42) 
(43) 

(44) 
(45) 



Using (|10l) , (|il|), (|39|) and ([13D , we can obtain expressions in the original 
coordinates. Taking the values of the above expressions at the disk and 
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using 



fl2l|) we compute the energy density and pressures 
"2 7 Va 2 -I- a l + (a 2 -l)7 2 r 2 



Pr 



Pr 



Pr 



Po 



a 1 + (a 2 — l)r 2 

1 + (a 2 - l) 7 2 r 2 ' 
1 + (a 2 - l)r 2 



v 7 !^ 2 [1 + (a 2 - l)r 2 f< 



(46) 

(47) 
(48) 



where po is a positive constant and < r < 1. 
From (HBT) we obtain the conditions 



Q 



> 1, 



ft 



to have e > at r = 0, and 



ft 2 /3 2 -2ft 2 /3+(ft 2 + l) < 0, 



(49) 



(50) 



in order to have e > at r = 1. We see that 7 and a can not be adjusted to 
have e > at the border of the disk and so e(r) have a zero at ro given by 



^0 



2a 2 ((3- 1) 



(3 2 -(1- « 2 ) (2/3-1) • 



(51) 



Now, to have e(r ) =0 with real values of 7 and a, we obtain the condi- 
tions 

1-r 2 _ „ 1 



a < 



1 ' 
r 2 



< < 



The value of j3 is given by 
1 



P 



a 2 rQ 



l + (a 2 -l)r 2 ±J(l-r 2 ) 2 -a 2 4 



(52) 



(53) 



Then, by choosing a value of r < l/\/2 we obtain values for a, (3 and 7 so 
that the solution satisfy the weak and strong energy conditions in the central 
region of the disk, for r < tq. 
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By choosing r 2 = 0.2 we can take a = 4, /3 — 5/4 and 7 = y5/3 so that 
the radial pressure is given by 



p r = p Vl - r 2 (l + 15r 2 ) 3 . (54) 

From (f£6|) - PS|) we can see that a = 2/3p r > so that the solution satisfy 
the weak and strong energy conditions in the central region of the disk. In 
Figure 2 we plot the functions e, p r , p v and a in the interval < r < 1 for 
this choice of the parameters, in p Q units. As in the first solution, the disk 
have a central region that satisfy the weak and strong energy conditions and 
a border region where e < 0, in violation of the weak energy condition. 



3.3 Bonnor and Sackfield Disks 

Finally, the third family of solutions is obtained by choosing < a < 1 in (|]) 
and transforming the disk into a oblate spheroidal shell. To solve the Laplace 
equation (|TJJ) we introduce oblate spheroidal coordinates ((, 77), related to the 
Weyl coordinates by 

^ 2 = ^ 2 (C 2 + l)(l-r7 2 ), Z = k( V , (55) 

where n = Vl — ct 2 , < £ < 00 and < 77 < 1. The shell is at £ = 
a/y/l - a 2 > 0. 

The general solution for the exterior of the shell, continuous across the 
disk, can be written as 

00 

$ = -]T2C 2 „P 2n (77)g 2ri (0, (56) 

n=0 

where (?2n(C) = ^ 2n+1 Q2n{.K)- The simplest case is given by 

$ = -27^-^, A = - 7 2 ln|^^, (57) 

where 7 is a real constant. As in the above two solutions, $ and A are 
continuous on the disk and vanish at infinity. The above solution was found 
by Zipoy |28| and Voorhees J2J| and interpreted by Bonnor and Sackfield [|TJ 
as the gravitational field of a pressureless disk. We shall see that the method 
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based in the conformal transformation not only leads to a nonzero radial 
pressure, but also to a nonzero azimuthal pressure. 

From the metric tensor in oblate spheroidal coordinates 



gtt 



9<p<f 



—exp{—2jcot 
R 2 exp {2y cot" 1 C} , 



\W 



C 2 + v 2 
l + C 2 



exp 1 27 cot 1 C| 



(58) 
(59) 

(60) 
(61) 



( |22|) - (gj) we compute the energy density and pressures 

"27VI - a 2 - a l + (l-a 2 )7 2 r 2 



we can obtain, using (|TUD, QTT]) , ( [55] ) and (JT5|), its expression in the original 
coordinates. Taking the values of the above expressions at the disk and using 



(62) 
(63) 
(64) 



Pr 



Pr 



Pr 



a 



(a* 



IV 2 



2^2 



1 + (1 - a z )Yr 
1 + (a 2 - IF 



p V 7 ! - r 2 



[l + (a 2 -l)r 2 ](« 2 + 1 )■ 



where po is a positive constant and < r < 1. 
From (|62|) we obtain the condition 



7\/i — a 2 



> 1, 



a 



(65) 



to have e > at r = 0. As in the second family of solutions, the condition 
(50) is imposed in order to have e > at r — 1. Again, 7 and a can not be 
adjusted to have e > at the border of the disk and so e(r) have a zero at r 
given by (|5l"l) . Now, to have e(r ) = with real values of 7 and a, we obtain 



a < 



- < ri < 1, 



(66) 
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so that the value of j3 is again given by (|5^). Then, by choosing a value of 
1 / \pl < r < 1 we obtain limits for 7 and a so that the solution satisfy the 
weak and strong energy conditions in the central region of the disk, r < r Q . 



By choosing tq 2 = 0.8 we can take a = 1/4, (3 — 5 and 7 = y 5/3 so that 
the radial pressure is given by 



Pr = Tj i^Ti- ( 67 ) 

From ( pf) - (f|8|) we can see that a = 2(3p r > and so again we have a solution 
that satisfy the weak and strong energy conditions for < r < r . In Figure 
3 we plot the functions e, p r , p v and a in the interval 0.85 < r < 1 for this 
choice of the parameters, in p Q units. As in the above two solutions, the 
disk have a central region that satisfy the weak and strong energy conditions 
and a border region where e < 0, in violation of the weak energy condition. 
However, in this case the central well behaved region can be extended to 
practically the entire disk by an appropriate selection of the parameters 7 
and a. 



4 Statics disks with Radial Tension 

We now find some solutions with radial tension by taking the negative sign 
of the square root in (||). There is only one family of solutions, corresponding 
to the case a = 1, obtained by transforming the disk into a spherical shell of 
radius one and the outer region of the disk into the inner region to the shell. 

To solve the Laplace equation (|TJ) we again introduce the spherical co- 
ordinates (|27D . The general solution for the Laplace equation ([14]) in the 
interior of the shell, continuous across the disk, can be written as 

00 

$ = -^2C 2n ft 2n P 2n (cos#), (68) 

n=0 

where C2 n are constants. The simplest case is obtained by taking the first 
term of the series and is given by 

$ = -2 7 , A = 0, (69) 

where 7 is a real constant. 
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The metric tensor is given by 

9tt = 



(70) 



g w = R 2 e 2 \ (71) 
g rr = \W'\ 2 e 2 \ (72) 

9rr = 9zz, (73) 



where R(r, z) and |W(a>)| are obtained from ( |10"D and (pTTf) . 

From ([TBI) and (|17D , we see that -R ab C£ 2 = and so we have a spacetime 



defect [pi. The energy density and tensions are given by 



e = - 2p v = -2p r = 2p Vl-r 2 , (74) 

where po is a positive constant and < r < 1. In this case we have a = 
and we can see that this solution satisfy all the energy conditions everywhere 
on the disk. This disk can be interpreted as being formed by cosmic strings 
located along the radial direction interwoven with loops. Another disklike 
topological defect, an infinite disk with zero radial stress, was presented in 
[ p~2]j . The energy density and the azimuthal pressure were given by e = 
p v — 1, also in agreement with the energy conditions. 

A second solution is obtained by taking the second term of the series in 
We find 



I = 1 (R 2 -2Z 2 ), A = ^-(R 2 -8Z 2 ), (75) 

where 7 is a real constant. As in the other solutions, <3> and A are continuous 
across the disk and vanish at infinity, where R — Z — 0. 
The metric tensor in Weyl coordinates is given by 

gtt = -exp{ 1 R 2 -2 1 Z 2 }, (76) 
g w = R 2 exp{2 1 Z 2 -jR 2 }, (77) 
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g rr = |^'| 2 exp |lA_4 7 2 jR 2 z 2_ 7jR 2 + 27Z 2| ) ^ 
9zz = 9rr- (79) 

Expressions in the original coordinates are obtained using (|T0|), (|TTD, and 

(PD- 

Taking the values of the above expressions at the disk and using ( |2"2"| ) - 
(|24|) we compute the energy density and pressures, 



~Pr 



P V = -Pr 



97V 4 - (6 + 87)7^ + 2(2 7 + 1)] , (80) 
87V - 97V - l] , (81) 



p r = —po exp I (3 + 47)7r 2 — 4.57 2 r 4 | Vl — r 2 , (82) 

where po is a positive constant and < r < 1. 

From (|80|) we can see that the weak energy condition is equivalent to 

97V 4 - (6 + 87)7r 2 + 2(2 7 + 1) > 0, (83) 

for < r < 1. This expression have two zeros at 

2 3 + 4 7 ± VI67 2 - 12 7 - 9 



" 9 7 



(84) 



and so there are three different classes of solutions corresponding to the values 
of 7 such that the above equation have either two, one or any real positive 
root. 

When 7 < —0.5 there is a solution r < 1 of ( |8"4| ) in such away that 
the disks have central regions with e < and border regions with e > 0. 
When -0.5 < 7 < 3(1 + y/E)/S we can see that e > everywhere on the 
disks, accordingly the weak energy condition is satisfied. Finally, when 7 > 
3(1 + v / 5)/8, there are two different real positive solutions of (H), r\ and 
r2, so that the disks have a central region with e > 0, when r < ri, a region 
with e < 0, when 77 < r < r2, and a region with e > 0, when r > T2- 

From ( P0| ) - (^) we can see that 

a = _ m .(4_6r 2 ). (85) 
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Hence <r = 0atr 2 = 2/3. Accordingly, the solutions are in agree with the 
strong energy condition only in a region, the central region of the disks when 
7 > 0, or the border region of the disks when 7 < 0. We also see from (PT|) 
that the azimuthal stress is zero at 

r± = ^ • (86) 

When I7I > 0.75 there are two different positive roots, r_ < r + < 1. Hence 
there is a region of the disks where the azimuthal stress is positive, for r_ < 
r < r + . On the other hand, when I7I < 0.75 the azimuthal stress is negative 
everywhere. 

In order to show the behavior of the different classes of solutions that can 
be obtained we consider some particular values of 7 . In Figure 4 - 7 we plot 
the functions e, p r , p v and a in the interval < r < 1 for 7 = —1, 7 = —0.5, 
7 = 1 and 7 = 2, respectively, in p units. We see that very different classes 
of solutions can be obtained, some of them with either the center of the disk 
or the border in accord with the energy conditions. 



5 Concluding Remarks 

In this work a method was presented which can be used to obtain static thin 
disks models with a well behaved nonzero radial stress. The general expre- 
ssion that was obtained for the radial stress admits solutions with positive 
stress, representing disks with radial pressure support, and solutions with 
negative stress, representing disks with radial tension support. Is easy to 
shown that, with the simple solution (|9]), when we consider solutions with 
positive stress is not possible to fulfill all the energy conditions at the entire 
disk, whereas that with negative stress we can obtain solutions that satisfy 
all the energy conditions everywhere on the disks, as the first solution studied 
at Sec. 4. 

The method that was described can be used to obtain many static disks 
solutions with a radial pressure support that presents a very reasonable phy- 
sical behavior and perhaps, by considering another more involved solutions of 
the Laplace equation (0), the solution can be adjusted to fulfill the weak and 
strong energy conditions everywhere on the disks, leading so to a physically 
acceptable source which would produce the considered vacuum solutions. 
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Also, as we will do in a future paper, the method can be extended to consider 
rotating thin disks, which have a stronger astrophysical interest. 
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